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Abstract 

The dependence of dijet production on the virtuality of the exchanged photon, 
Q 2 , has been studied by measuring dijet cross sections in the range ~ Q 2 < 
2000 GeV 2 with the ZEUS detector at HERA using an integrated luminosity 
of 38.6 pb" 1 . Dijet cross sections were measured for jets with transverse energy 
Ej, > 7.5 and 6.5 GeV and pseudorapidities in the photon-proton centre-of-mass 
frame in the range —3 < if et < 0. The variable x° bs , a measure of the photon 
momentum entering the hard process, was used to enhance the sensitivity of 
the measurement to the photon structure. The Q 2 dependence of the ratio of 
low- to high-x° bs events was measured. Next-to-leading-order QCD predictions 
were found to generally underestimate the low-x° bs contribution relative to that 
at high £° bs . Monte Carlo models based on leading- logarithmic parton-showers, 
using a partonic structure for the photon which falls smoothly with increasing 
Q 2 , provide a qualitative description of the data. 
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1 Introduction 



Interactions involving real or quasi- real photons (Q 2 ~ 0, where Q 2 is the virtuality 
of the photon) are well described by calculations that use a partonic structure for the 
photon [1]. However, in deep inelastic scattering (DIS), where Q 2 is large, the virtual 
photon is commonly treated as a point-like object and used as a probe of the partonic 
structure of nucleons [2]. In this paper, dijet production is investigated over a large range 
of incident photon virtualities, including photoproduction, DIS, and the transition region 
between them. Both the HI [3] and ZEUS [4] collaborations have previously studied the 
transition between photoproduction and DIS by measuring inclusive jet and dijet cross 
sections in ep collisions. 

Two processes contribute to the jet photoproduction cross section at leading order (LO) 
in quantum chromodynamics (QCD) [5,6]: direct, in which the photon couples as a point- 
like particle to quarks in the hard scatter; and resolved, in which the photon acts as a 
source of partons. Both processes can lead to two jets in the final state. The x° bs variable, 
which is the fraction of the photon momentum participating in the production of the dijet 
system, is used to separate the two processes since resolved (direct) processes dominate 
at low (high) x° bs values [7]. 

In conventional fixed-order QCD calculations, only point-like photon interactions con- 
tribute to jet production in DIS. However, two scales play a role in the interaction: Q 
and the jet transverse energy, £f\ For high Q 2 {Q 2 > (E^) 2 ), QCD predicts that the 
photon will behave as a point-like object. For Q 2 <C (Eji 1 ) 2 , the photon may have an 
effective partonic structure, even for relatively large values of Q 2 , which is resolved at 
a scale related to the transverse energy of the jets. Therefore, resolved processes may 
contribute significantly to the jet cross section. The ratio of cross sections evaluated in 
different x° bs ranges is particularly sensitive to the resolved component. 

In this paper, the validity of the above approaches in photoproduction and DIS is studied 
by measuring dijet cross sections differential with respect to Q 2 , E^ tl and rj F , where 
E^ a is the Ex of the jet in the accepted rapidity range which has the highest transverse 
energy, and rf is the pseudorapidity of the most forward jet. The ratio of low- to high-£° bs 
components is measured as a function of Q 2 in different regions of E T , where E T is the 
square of the average transverse energy of the two jets with highest transverse energy. 

The data used in this analysis correspond to an integrated luminosity six times larger 
than that used in the previous ZEUS study [4]. Next-to- leading-order (NLO) QCD cal- 
culations [8-11] have been compared to measurements that span a large range of photon 
virtualities. The predictions of leading-logarithm parton-shower (PS) Monte Carlo (MC) 
models are compared to the data in the transition region between photoproduction and 
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DIS, where current NLO calculations are not applicable. 



2 Experimental set-up 

The data were collected during the 1996 and 1997 running periods, when HERA operated 
with protons of energy E p = 820 GeV and positrons of energy E e = 27.5 GeV, and 
correspond to an integrated luminosity of 38.6 ± 0.6 ptr 1 . 

The ZEUS detector is described in detail elsewhere [12]. The most important components 
used in the current analysis were the central tracking detector (CTD), the uranium- 
scintillator calorimeter (CAL) and the beam pipe calorimeter (BPC). 

Charged particles are tracked in the central tracking detector (CTD) [13], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 
polar-angle 1 region 15° < < 164°. The transverse-momentum resolution for full-length 
tracks is a(p T )/p T = 0.0058p T © 0.0065 © 0.0014/p r , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [14] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18/ \[E for electrons 
and a(E)/E = 0.35/ \[E for hadrons, with E in GeV. 

The BPC [15] was installed 294 cm from the interaction point in the positron direction 
in order to tag scattered positrons at small angles with respect to the positron beam 
direction (15-34 mrad). It measured both the energy and impact position of the scattered 
positron at the BPC surface. The relative energy resolution of the BPC is 0.17/v^E and 
the position resolution is 0.5 mm. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > e'yp, 
where the photon was measured with a lead-scintillator calorimeter [16] at Z = —107 m. 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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3 Theoretical framework 



In photoproduction, perturbative QCD (pQCD) calculations of dijet cross sections can 
be written as a convolution of the subprocess cross section with the parton distribution 
functions (PDFs) of the photon and proton: 

1 1 

^°ep->ejctjct = J dyf^/ e (y,fj,ji) J dx~ / f a /~ / (x.y,fi R ,fi Fl ) x 

a,b 

1 



where y, x 7 and x v are the longitudinal momentum fractions of the almost-real photon 
emitted by the positron, the parton a in the photon and the parton b in the proton, 
respectively. The function / 7/e is the flux of photons from the positron, and f a / y (/&/ p ) 
represents the PDF of parton a (b) in the photon (proton). The factorisation scale for the 
photon (proton) is denoted by /xp 7 (lif p ) an d llr represents the renormalisation scale. The 
subprocess cross section, da a ^j et j et , describes the short-distance structure of the interac- 
tion. For direct processes in the above formula a is replaced by 7 and / a / 7 (a; 7 , Q 2 , /if> 7 ) is 
given by 5(1 — x 7 ). 

In DIS, the photons are virtual (7*) and usually only direct processes are considered. 
Effective resolved terms appear only as higher-order corrections. 

In the transition region between DIS and photoproduction, a virtual-photon structure [17- 
19] may be introduced. In general, the virtual-photon PDFs / a / 7 * contain two terms, 

/a/7* (^7*) Q 2 i A^*) = f a / ? ( X 7* ' Q 2 1 A*F 7 * ) + /a/ 7 * ( X 7* ' Q 2 ' ) ' 

the first associated with the non-perturbative hadronic component (y non -P crt ) ; j n which 
the photon fluctuates into an intermediate meson-like hadronic state, and the second / pert , 
unique to the photon, which expresses the coupling of the photon to a high-virtuality qq 
pair, calculable in pQCD. Perturbative QCD predicts that the contribution to the dijet 
cross section from resolved processes should decrease relative to the contribution from di- 
rect processes as the virtuality of the photon increases towards [xr. The non-perturbative 
component of the virtual-photon PDFs decreases as Q~ 4 , whereas the perturbative com- 
ponent decreases as \b.(li r /Q 2 ). 

Two parameterisations of the virtual-photon PDFs, SaS [17] and GRS [19], are available. 
Both are extrapolations of the real-photon PDFs to the virtual-photon regime. They differ 
in the treatment of the non-perturbative component. In the case of the SaS sets, a fit to 
a coherent sum of the lowest-lying vector-meson states p, u and <fi has been performed, 
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whereas, in the case of GRS, the non-perturbative part has been estimated using the 
PDFs of the pion. 



4 Cross section definition 

Dijet cross sections differential in Q 2 , E^ 1 and r] F were measured. The ratios of cross 
sections for low (< 0.75) to high (> 0.75) x° bs are presented. The variable x° bs is defined 
as 

obs _ ^jetsV 



7 Ehadrons(^ ~ Vz) 

where E' iet and are the energy and the longitudinal momentum of the jet. The upper 
sum runs over the two jets with highest transverse energy and the lower sum runs over 
all final state hadrons. 

The cross sections were measured in the range ~ Q 2 < 2000 GeV 2 and 0.2 < y < 0.55. 
Jets were reconstructed with the kr cluster algorithm [20] applied in the photon-proton 
centre-of-mass frame, in the longitudinally invariant inclusive mode [21]. At least two jets 
were required within the pseudorapidity range —3 < rf ct < 0, satisfying E^ 1 > 7.5 GeV 
and E% t2 > 6.5 GeV. 



5 Data selection and jet search 

A three-level trigger was used to select events online [12,22]. In the third-level trigger the 
events were required to have at least two jets with a transverse energy of E^ > 4 GeV 
and a pseudorapidity of rf et < 2.5 in the laboratory frame. 

The sample was separated offline into subsamples corresponding to three different Q 2 
ranges: 

• DIS sample: events were selected by requiring that the outgoing positron was measured 
in the CAL [23]. The energy of the scattered positron, E e >, was required to be above 
10 GeV, with 1.5 < Q 2 < 2000 GeV 2 ; 

• BPC sample: events at low Q 2 were selected by requiring that the scattered positron 
was measured in the BPC. These events were required to have E e > > 12.5 GeV and 
0.1 < Q 2 < 0.55 GeV 2 ; 

• Photoproduction sample: events were selected by requiring that the scattered positron 
was not observed in the CAL, implying Q 2 < 1 GeV 2 with a median Q 2 ~ 10~ 3 GeV 2 . 
A small fraction of this sample (0.6 %) is also contained in the BPC sample. 
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For all three samples, hadronic kinematic variables and jets were reconstructed using 
a combination of track and CAL information which optimises the resolution [24]. The 
selected tracks and CAL clusters are referred to as Energy Flow Objects (EFOs). 

The method reported in a previous publication [25] was used to correct the EFOs for 
energy losses in inactive material in front of the CAL. The jet-energy-scale uncertainty is 
within ±1% for E^ > 7.5 GeV and increases to ±3% for lower E^ values. 

Additional cuts, similar to those used in an earlier analysis [4], were applied offline to all 
samples: 

• a reconstructed event vertex consistent with the nominal interaction position was 
required, \Z vtx \ < 40 cm; 

• to suppress the background from events with a misidentified positron, the variable 
y e — 1 — 2E ^fj' cos9 ,) was required to satisfy y e < 0.8, where 9 e > is the polar angle of 
the scattered positron; 

• for the DIS sample, a fiducial volume cut was applied to the positron position (|V e | > 
14 cm or \Y e \ > 9 cm, where X e and Y e are the impact positions of the positron on 
the face of the CAL) in order to avoid the low-acceptance region adjacent to the rear 
beam pipe; 

• for the BPC sample, the reconstructed impact position on the BPC surface was con- 
strained to be within the fiducial- region of the BPC [15]; 

• for the photoproduction sample, events with a scattered-positron candidate in the 
CAL were rejected, as in a previous ZEUS analysis [26]; 

• all samples were required to satisfy 0.2 < ?/jb < 0.55, where yjB = J2i(Ei~Ezi)/2E e [27] 
is an estimator of y. The sum runs over all EFOs. Ezi = E^ cos 6i, where Ei is the 
energy of EFO i with polar angle d\ with respect to the measured Z-vertex of the 
event. The lower cut removes beam-gas events and the upper cut is imposed due to 
the restricted acceptance of the BPC detector. 

Prior to jet finding, the EFOs were boosted to the photon-proton centre-of-mass frame. 
In the DIS and BPC samples the boost was calculated using the reconstructed momentum 
of the scattered positron. In the photoproduction sample ?/jb was used in performing the 
boost. 

The Jzt cluster algorithm was applied to the boosted EFOs in the photon-proton centre- 
of-mass frame to reconstruct jets. At least two jets were required in each event within the 
pseudorapidity range — 3 < rf ct < and were ordered according to decreasing E^ . They 
were further required to satisfy E^ n > 7.5 GeV and E^ t2 > 6.5 GeV. After all cuts, the 
photoproduction/BPC/DIS sample contained 419911/2481/45100 dijet events. The BPC 
sample is a subset of the photoproduction sample. 
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6 Acceptance corrections 



The programs Herwig 5.9 [28] and Pythia 6.1 [29] were used to generate events for 
resolved and direct processes over the whole Q 2 range. Events were generated using 
GRV-LO [30] for the photon PDFs and MRSA [31] for the proton PDFs. To study 
the dependence of the acceptance corrections on the choice of photon and proton PDFs, 
the GRS-LO and CTEQ5M1 [32] parameterisations were used, respectively. In both 
generators, the partonic processes are simulated using LO matrix elements, with the 
inclusion of initial- and final-state parton showers. Hadronisation is performed using a 
cluster model [33] in the case of Herwig and the Lund string model [34] in the case 
of Pythia. For the measurements presented in this paper, the Herwig and Pythia 
programs were used to correct the data for acceptance. The corrections provided by 
Herwig were used as default values and those given by Pythia were used to estimate 
the systematic uncertainties associated with the treatment of the parton shower and 
hadronisation. 

All generated events were passed through the ZEUS detector and trigger simulation pro- 
grams based on Geant 3.13 [35]. They were reconstructed and analysed by the same 
program chain as the data. The jet search was performed using EFOs in the same way 
as for the data. The same jet algorithm was also applied to the final-state particles. The 
jets found in this way are referred to as hadronic jets. 

The acceptance corrections take into account the efficiency of the trigger, the selection 
criteria and the purity and efficiency of the jet reconstruction. The differential dijet cross 
sections were obtained by applying bin-by-bin corrections to the measured distributions. 
The predictions of the generators Herwig and Pythia for the uncorrected distributions 
were compared to the data for the above parameterisations of the photon and proton 
PDFs. The contributions from direct and resolved processes were added according to a 
fit to the uncorrected x° bs distribution in the data. A good description of the , rf et , 
Q 2 and y data distributions was given by both Herwig and Pythia. 

For the photoproduction sample the bin-by-bin correction factor was approximately 1.2. 
This increased to approximately 6 for the BPC sample due to the geometric acceptance of 
the BPC detector [15]. For 1.5 < Q 2 < 4.5 GeV 2 , the correction factor was approximately 
3 due to the fiducial volume cut (see section EJ). For Q 2 > 4.5 GeV 2 , the bin-by-bin 
correction factors differed from unity by less than 10%. 
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7 QCD calculations 



7.1 NLO calculations 

The NLO QCD calculations of jet production cross sections in DIS used in this analysis are 
based on the programs DlSASTER-| — h [8] and DlSENT [9]. In these programs, the photon 
is treated as a point-like probe. Contributions from hadron-like resolved processes are 
not included. They use the subtraction method [36] and the massless MS renormalisation 
and factorisation schemes. Their predictions agree to within ±3%. In Section El only the 
calculations using DisasterH — h are compared to the data because this program allows 
a wider parameter selection than DlSENT. In the calculations, the number of flavours 
was set to five. The renormalisation and factorisation scales were set to /i 2 = fi 2 R = 
I^f = Q 2 + (E^ 1 ) 2 or Q 2 , and a s (fif>) was calculated at two loops using = 226 MeV 
corresponding to as[Mz) — 0.118. The CTEQ5M1 sets were used for the proton PDFs. 

Many calculations of jet photoproduction at NLO exist [10,11,37-41], all of which agree 
to within (5 — 10)% [41,42]. The calculations of Frixione and Ridolfi [10,11] uses the 
subtraction method. In this calculation the number of flavours was set to five and the 
factorisation and renormalisation scales to ji 2 = (E^ 1 ) 2 . For the calculation of ck s (//r), 

(5) 

A^i = 226 MeV was used. For the proton PDFs, the CTEQ5M1 sets were used, and for 
the real photon PDFs the GRV and AFG [43] parameterisations were used. 

Samples of events generated using the Heracles 4.6.1 [44] MC program with the Djan- 
GOH 1.1 [45] interface to the hadronisation programs were used to estimate hadroni- 
sation corrections for the NLO QCD predictions calculated using Disent and Disas- 
terH — h The QCD cascade is simulated using the colour-dipole model [46] including 
the LO QCD diagrams as implemented in Ariadne 4.08 [47] or with the MEPS model 
of Lepto 6.5 [48]. Both Ariadne and Lepto use the Lund string model [34] for the 
hadronisation. For the photoproduction NLO prediction, the Herwig and Pythia MCs 
were used to estimate the hadronisation corrections. 

First-order QED radiative effects were also estimated using Heracles and found to be 
1% or less. Corrections for these effects have not been applied to the NLO calculations. 

The predictions to be compared with the data were corrected for hadronisation effects 
using a bin-by-bin procedure according to da = da NLO ■ C^ d , where <i<r NLO is the cross 
section for partons in the final state of the NLO calculation. The hadronisation correction 
factor was defined as the ratio of the dijet cross sections before and after the hadronisation 
process, Chad = ^ a Mc° n V^ a Mc r ° ns - The value of Chad was taken as the mean of the ratio 
obtained using the predictions of two different generators (Ariadne and Lepto for DIS, 
and Herwig and Pythia for photoproduction) and was found to lie between 1.1 (large 
Q 2 ) and 1.2 (small Q 2 and photoproduction). 
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7.2 Monte Carlo predictions 

Predictions of Herwig 6.4 [49] using CTEQ5L for the proton PDFs and SaS2D for the 
photon PDFs were generated using parameters tuned [50] to many previous HERA and 
LEP measurements. In the SaS2D parameterisation the structure of the virtual photon 
is suppressed with increasing Q 2 . Predictions were also generated with this suppression 
switched off. 

8 Systematic uncertainties 

8.1 Experimental uncertainties 

A detailed study of the sources contributing to the systematic uncertainties of the mea- 
surements was performed. This study includes (a typical contribution to the uncertainty 
in the cross section for each item is indicated in parentheses): 

• using the Pythia generator to evaluate the acceptance corrections to the observed 
dijet distributions (+6%); 

• using different parameterisations of the photon (GRV-LO and GRS) and proton (MRSA 
and CTEQ5M1) PDFs for the generation of the Herwig MC samples (±2%); 

• varying the Ej, cut by the resolution (±8%); 

• varying the other selection cuts by their respective resolution (< ±2%) 

• adding the contributions from direct and resolved processes according to the default 
cross sections as predicted by Herwig (—3%). 

All the above systematic uncertainties were added in quadrature. The effect of the uncer- 
tainty in the absolute energy scale of the jets on the dijet cross sections was approximately 
±9% at low Q 2 , decreasing to ±6% at high Q 2 . This uncertainty is highly correlated and is 
shown separately in the figures. In addition, there is an overall normalisation uncertainty 
from the luminosity determination of 1.6%, which is not shown in the figures. 

8.2 Theoretical uncertainties 

The NLO QCD predictions for the dijet cross sections are affected by the following theo- 
retical uncertainties (typical values for the uncertainties are quoted): 

• uncertainties due to terms beyond NLO, estimated by varying \i by factors 2 and 0.5 
(20% at low Q 2 and 7% at high Q 2 , in the case of /x 2 = Q 2 + (E^) 2 ); 
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• uncertainties in the hadronisation corrections, estimated as half the spread between 
Chad values obtained using the Herwig, Pythia, Lepto and Ariadne models (2 - 

3%); 

• uncertainties on the calculations due to as and the proton PDFs, estimated by using 
the MRST sets of parameterisations (5%). These uncertainties were cross-checked 
using an alternative method [51], which uses the covariance matrix of the fitted PDF 
parameters and derivatives as a function of x. 

The above theoretical uncertainties were added in quadrature to give the total uncertainty 
on the predictions. 

9 Results 

9.1 Single-differential dijet cross sections 

Figure H shows the differential dijet cross section, da/dQ 2 , for > 7.5 GeV, E^ t2 > 
6.5 GeV, -3 < r/ jct < 0, 0.2 < y < 0.55 and 0.1 < Q 2 < 2000 GeV 2 together with 
the photoproduction measurement. The cross section split in the direct-enhanced region 
(x° bs > 0.75) and the resolved-enhanced region (x° bs < 0.75) is also shown. All the cross 
sections are given in the tables ^ 121 and El 

The measurements cover a wide range in Q 2 , including the transition region from pho- 
toproduction to DIS. The measured cross sections fall by about five orders of magni- 
tude over this Q 2 range. The cross section for x° bs < 0.75 falls more rapidly than that 
for x° bs > 0.75. Even though the dijet cross section is dominated by interactions with 
x° bs > 0.75 for Q 2 ~ 10 GeV 2 , there is a contribution of approximately 24% from low-:c° bs 
events with Q 2 ~ 500 GeV 2 . 

The NLO QCD calculations are compared to the measured da/dQ 2 in Figs, d and 121 
The prediction 2 with /i 2 = Q 2 + (E^ 1 ) 2 , shown in Fig. ^ describes the shape of the 
measured total dijet cross section but underestimates its magnitude by approximately 
30%. The renormalisation scale uncertainty was evaluated also for the low- and high-x° bs 
cross sections. For the high-x° bs cross section this uncertainty was similar to that on the 
total cross section. In the case of the low-x° bs cross section, the uncertainty was almost 
constant at around ±30%. Taking these uncertainties into account, the measured cross 
section for x° bs > 0.75 is reasonably well described by the calculation shown in Fig. for 
all Q 2 . However, the prediction dramatically underestimates the measured cross section 
for xf < 0.75. 

2 The two lowest Q 2 bins are outside the range of applicability of the DisasterH — h program. 
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The prediction with /i 2 = Q 2 is shown in Fig. |21 It has a much larger renormalisation-scale 
uncertainty than the prediction using /i 2 = Q 2 + (E^) 2 , and within this uncertainty it is 
consistent with the data. 

A possible explanation of the disagreement, and for the large uncertainties in the pre- 
diction at low-x° bs values, is that effects arising from the structure of the photon are 
expected in this region, whereas the contribution predicted by DisasterH — h comes only 
from large-angle particle-emission diagrams included in the NLO corrections to the dijet 
cross section. 

In photoproduction, the low-x° bs component of the data becomes dominant. The photo- 
production measurement is well described by the photoproduction NLO prediction, using 
the GRV photon PDF. 

9.2 Double-differential dijet cross sections 

The dijet cross section, d 2 a/dQ 2 dE^ tl , as a function of E^ 11 in different Q 2 ranges is shown 
in Fig. El and given in the Tables 0] and The measurements extend up to transverse 
energies of approximately 40 GeV. The E^ tl distribution falls less steeply as Q 2 increases. 
Figure El also shows the NLO QCD predictions. The NLO calculation for photoproduction 
using GRV for the photon PDFs gives a good description of the E^ 1 cross section. At 
higher Q 2 , the calculation using /i 2 = Q 2 + (E^ 1 ) 2 is in agreement with the data for 
the lowest and highest jet transverse energies, but lies below the data for intermediate 
E^ 1 values. The prediction with fi 2 = Q 2 again agrees with the data, within the large 
theoretical uncertainties (not shown). 

The differential cross-section d 2 a/dQ 2 dr] F as a function of if is shown in Fig. 0] for 
different ranges of Q 2 and given in the Tables |H1 and The cross section as a function of 
7] F is more sensitive to the resolved photon component in the forward direction 3 . In all 
Q 2 regions, the measured cross section increases with rf in the region —2.5 < r] F < —1.5. 
For 7] F > —1.5, the cross section decreases as r] F increases for Q 2 ~ 10 GeV 2 , whereas 
in photoproduction and at low Q 2 the cross section increases. The NLO prediction for 
photoproduction describes the measured cross section. At low Q 2 the NLO prediction 
using jj? = Q 2 + (E^ t ) 2 underestimates the measured cross section in the forward direction. 
The prediction with fi 2 = Q 2 again agrees reasonably well with the data within large 
theoretical uncertainties (not shown). 

The differences between the data and NLO calculations may be due to the persistence of 
a resolved component at Q 2 > 1 GeV 2 . To study this in more detail, the ratio of dijet 

3 Since r\ here is defined in the hadronic centre-of-mass frame, the forward region in the laboratory frame 
corresponds to r\ > — 1. 
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cross sections for high and low x° bs values is presented in the next subsection. 
9.3 Ratios of dijet cross sections 

The Q 2 dependence of the direct- and resolved-enhanced components of the dijet cross 
section has been studied in more detail using the ratio 

a(xf s < 0.75) 
~~ cr(z° bs > 0.75) ' 

A number of experimental and theoretical uncertainties cancel in this ratio, so that the 
presence of a resolved contribution can be investigated at higher precision than in the 
individual cross sections. 

Figures 03 and El show the ratio R as a function of Q 2 in three different regions of E T . 

2 2 

The Q 2 dependence of the data is stronger at low E T than at higher E T , showing that 

2 2 

the low-x° bs component is suppressed at low Q 2 as E T increases, and at low E T as Q 2 
increases. The ratio is also given in Table |H1 

Predictions of the Herwig MC program using the SaS2D parameterisation of the photon 
PDFs are compared to the data in Fig. E3 The SaS2D parameterisation contains the 
suppression of the virtual photon structure with increasing Q 2 . The predictions fall with 
increasing Q 2 and qualitatively reproduce the data. However, the predictions using SaS2D 
with the suppression of the virtual photon structure switched off are relatively constant 
with Q 2 . 

The NLO calculations are compared to the data in Fig. The photoproduction calcula- 
tions using GRV are in reasonable agreement with the data, whereas those using AFG are 

2 

below the data. In the DIS region, the predictions lie below the data at low E T . However, 
some suppression in the ratio as a function of Q 2 is observed. 



10 Summary and conclusions 

Dijet differential cross sections have been measured in the range £ Q 2 < 2000 GeV 2 
with 0.2 < y < 0.55, -3 < r/ jct < and £^ tl,jct2 > 7.5 and 6.5 GeV, as a function of Q 2 , 
E^ a and r] F in the photon-proton centre-of-mass frame. These precise measurements, 
spanning a large range of photon virtualities, including photoproduction, DIS, and the 
transition region between them, are qualitatively described by leading-logarithmic parton- 
shower MC models which introduce virtual photon structure, suppressed with increasing 
Q 2 . These data may constrain such parton densities significantly if used in future fits. 
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The currently available next-to-leading-order QCD calculations have large uncertainties at 
low Q 2 , where the presence of a resolved-photon contribution may be expected. Improved 
higher-order or resummed calculations are needed. In DIS, the NLO QCD predictions 
generally underestimate the cross section at low-x° bs relative to that at high x° bs . 
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Q bin (GeV ) 


da/dQ' A stat A syst A ES (pb/GeV^) 


0, 1 


_ _ _ „ , . - _ _i_i no i o i v 

9280 ±113 ±?JSo 


0.1, 0.55 


2250 ±45.2 —2io -227 


1.5, 4.5 


167 ±2 22 + 14 ' 1 +12 - 5 
10/ zrz.zz _ lig _ 152 


4.5, 10.5 


54.5 ±0.54 ^ +22 


10.5, 49 


11.9 ±0.093 t°; 3 l ±?S 


49, 120 


2.27 ±0.027 1°;^ 


120, 2000 


nnqc; +0 001 1 +0-0019 +0.0059 
U.UJO ICU.UUII _o.0048 -0.0068 



Table 1: Measured dijet cross-sections da/dQ 2 . The statistical, systematic and 
jet energy scale, Aes, uncertainties are shown separately. 



Q 2 bin (GeV 2 ) 


da/dQ 2 A stat A syst (pb/GeV 2 ) 


0, 1 


5710 ±96.4 ttl 7 


0.1, 0.55 


1270 ±32.3 tux 


1.5, 4.5 


87.8 ±1.54 till 


4.5, 10.5 


24.5 ±0.35 tlH 


10.5, 49 


4.21 ±0.051 t° 5 27 


49, 120 


0.72 ±0.015 t° i 7 


120, 2000 


0.022 ±0.00052 tomll 



Table 2: Measured dijet cross-sections da/dQ 2 for x° bs < 0.75. The statistical 
and systematic uncertainties are shown separately. 



Q 2 bin (GeV 2 ) 


da/dQ 2 A stat A syst (pb/GeV 2 ) 


0, 1 


3620 ±65.8 tlti 


0.1, 0.55 


980 ±31.9 tlri 


1.5, 4.5 


79.5 ±1.61 ttf 6 


4.5, 10.5 


30 ±0.42 t 2 °i 


10.5, 49 


7.74 ±0.08 t i 7 


49, 120 


1.63 ±0.024 toll 


120, 2000 


0.077 ±0.0011 ±8:oo5? 



Table 3: Measured dijet cross-sections da/dQ 2 /or x° bs > 0.75. The statistical 
and systematic uncertainties are shown separately. 
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E% n bin (GeV) 


d 2 a/dE^dQ 2 A stat A syst A ES (pb/GeV 3 ) 


< Q 2 < 1 GeV 2 


7.5, 10 


1740 ±34.8 tlTs ii£ 


10, 13 


1010 ±20.4 tlH -lis 


13, 17 


007 10 c +21.5 +35.7 
00 ' =no.-J _g. 51 _39 9 


17, 22 


76 4 ±3 22 +3 ' 9 + 9 - 45 


22, 29 


ic + 1 9 c +0.1 +1.96 
ld 3=1. ZO -0.085 -2.03 


29, 50 


1.75 ±0.092 t^ 7 l 9 ; 22 


0.1 < Q 2 < 0.55 GeV 2 


7.5, 10 


484 ±22.6 i^; 3 1^ 


10, 13 


243 ±12.1 ±^ _}f 


13, 17 


68.7 ±4.91 t1| 2 ±£ 22 


17, 22 


17.9 ±2.19 ±}; 33 ±842 


22, 29 


2.49 ±0.67 ±°;1 l^ 3 


29, 50 


n 00 _un 1 7 +0.071 +0.048 
0.38 ±0.17 _ 083 _ 


1.5 < Q 2 < 4.5 GeV 2 


7.5, 10 


35.2 ±0.72 tl'rt 


10, 13 


17.6 ±0.41 I / 96 1?^ 


13, 17 


cnq xnifi +0.52 +0.36 
O.ua H=U.1U -0.044 -0.26 


17, 22 


1.64 ±0.091 t° if T H 


22, 29 


n -4-n n<« +0.0087 +0.015 

U.OZ 3ZU.UOO -0.035 -0.01 


29, 50 


n +n nnfn + 0023 + - 0024 

U.UCH ICU.UUUO -0.0049 -0.0035 


4.5 < Q 2 < 10.5 GeV 2 


7.5, 10 


10.2 ±0.16 I?; 23 ±l 21 


10, 13 


5.86 ±0.1 +_° il 


13, 17 


2.08 ±0.048 ±g:8S ±S:S 


17, 22 


n ^9 -4-n no +0.0084 +0.034 
u - dZ mu.uz -0.033 -0.019 


22, 29 


0.13 ±0.0097 t° Zl 


29, 50 


01 1 ±0 001 5 +0 - 002 +0-00039 

U.U11 ICU.UUIO -0.00084 -0.00059 



Table 4: Measured dijet cross-section d 2 a/dQ 2 dE^ tl . 
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E^ 1 bin (GeV) 


d 2 a/dE^dQ 2 A stat 


Asyst 


A ES (pb/GeV 3 ) 


10.5 < Q 2 < 49 GeV 2 


7.5, 10 


1.9 


±0.024 


+0.15 
-0.059 


+0.17 
-0.25 


10, 13 


1.36 


±0.018 


+0.057 
-0.041 


+0.091 
-0.082 


13, 17 


0.51 


±0.0091 


+0.024 
-0.013 


+0.028 
-0.025 


17, 22 


0.14 


±0.0041 


+0.0029 
-0.007 


+0.0031 
-0.0086 


22, 29 


0.036 


±0.0018 


+0.0014 
-0.0021 


+0.0017 
-0.0012 


29, 50 


0.0031 


±0.00031 


+0.00018 
-0.00011 


+0.00028 
-0.00034 


49 < Q 2 < 120 GeV 2 


7.5, 10 


0.34 


±0.0071 


+0.025 
-0.024 


+0.038 
-0.053 


10, 13 


0.23 


±0.0049 


+0.014 
-0.011 


+0.018 
-0.016 


13, 17 


0.11 


±0.0031 


+0.0016 
-0.0071 


+0.0058 
-0.0047 


17, 22 


0.039 


±0.0016 


+0.00087 
-0.001 


+0.0011 
-0.0014 


22, 29 


0.0059 


±0.00047 


+0.001 
-0.00023 


+0.00031 
-0.00018 


29, 50 


0.00072 


±8.5 x 10" 5 


+0.00023 
-3xl0" 5 


+8.3X10" 5 
-5.2xl0- 5 


120 < Q 2 < 2000 GeV 2 


7.5, 10 


0.01 


±0.00025 


+0.00041 
-0.00098 


+0.001 
-0.0014 


10, 13 


0.0099 


±0.00021 


+0.00043 
-0.00076 


+0.00059 
-0.0006 


13, 17 


0.0056 


±0.00014 


+0.00011 
-0.00023 


+0.00022 
-0.00025 


17, 22 


0.0021 


±7.7 x 10" 5 


+6.7X10- 5 
-0.00011 


+0.00014 
-6.4xl0- 5 


22, 29 


0.00069 


±3.6 x 10~ 5 


+7.1xl0~ 5 
-3.9xl0- 5 


+1.8X10 -5 
-3.3xl0" 5 


29, 50 


6 x 10~ 5 


±5.2 x 10~ 6 


+5xl0" 6 
-4.2xl0- 6 


+3.4xl0~ 6 
-2xl0- 6 



Table 5: Measured dijet cross-s ection d 2 ojdQ 2 dEf x . 
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rj F bin 


d 2 a/d v F dQ 2 A stat A syst A ES (pb/GeV 2 ) 


< Q 2 < 1 GeV 2 


-3, -1.8 


771 +97 r +39.1 +92.8 
111 =cz/ -38.3 -119 


-1.8, -1.4 


3790 ±108 IJj" till 


-1.4, -0.8 


4510 ±96.9 t^ 9 tf 7 ? 


-0.8, 


conn _i_qq +66.6 +494 

dOUU ICJJ.O _229 -606 


0.1 < Q 2 < 0.55 GeV 2 


-3, -1.8 


229 ±18.1 tlH tlli 


-1.8, -1.4 


939 ±64.2 tlli ttli 


-1.4, -0.8 


1070 ±55 t™ tlH 


-0.8, 


1240 ±58.7 l|i tZ 


1.5 < Q 2 < 4.5 GeV 2 


-3, -1.8 


15.3 ±0.62 1 2 q 8 ^ 2 Q5 


-1.8, -1.4 


75.8 ±2.47 ttil ±|;f 


-1.4, -0.8 


81.5 ±1.96 ±1° 15 ±|;^ 


-0.8, 


87.1 ±1.78 tH t 6 7 H 


4.5 < Q 2 < 10.5 GeV 2 


-3, -1.8 


6.67 ±0.19 ±g^ t° 7 sl 


-1.8, -1.4 


26 ±0.61 ±BiI ±5:S 


-1.4, -0.8 


26.8 ±0.47 tlil ±i:S 


-0.8, 


24.9 ±0.41 ^ l^ 2 



Table 6: Measured dijet cross-section d 2 a/dQ 2 dr] F . 
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t] F bin 


d 2 a/dr] F dQ 


2 A 


^syst 


A ES (pb/GeV 2 ) 


10.5 < Q 2 < 49 GeV 2 


-3, -1.8 


1.45 


±0.031 


+0.17 
-0.036 


+0.15 
-0.2 


-1.8, -1.4 


6.57 


±0.11 


+0.48 
-0.25 


+0.44 
-0.47 


-1.4, -0.8 


6.21 


±0.084 


+0.37 
-0.16 


+0.36 
-0.45 


-0.8, 


4.78 


±0.064 


+0.4 
-0.22 


+0.28 
-0.37 


49 < Q 2 < 120 GeV 2 


-3, -1.8 


0.3 


±0.01 


+0.032 
-0.042 


+0.036 
-0.045 


-1.8, -1.4 


1.27 


±0.034 


+0.14 
-0.08 


+0.098 
-0.1 


-1.4, -0.8 


1.24 


±0.026 


+0.074 
-0.08 


+0.076 
-0.09 


-0.8, 


0.87 


±0.018 


+0.065 
-0.056 


+0.059 
-0.065 


120 < Q 2 < 2000 GeV 2 


-3, -1.8 


0.0097 


±0.00035 


+9.2X10" 5 
-0.0011 


+0.00093 
-0.0014 


-1.8, -1.4 


0.055 


±0.0015 


+0.0011 
-0.0067 


+0.0035 
-0.0032 


-1.4, -0.8 


0.051 


±0.001 


+0.003 
-0.0012 


+0.0027 
-0.0032 


-0.8, 


0.039 


±0.00082 


+0.002 
-0.0019 


+0.0022 
-0.0024 



Table 7: Measured dijet cross-section d 2 o~/dQ 2 dr] F . 
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Q 2 bin 


R A sta t A syst A E S 


49 < E l T < 85 GeV 2 


0, 1 


919 in C\7K +0.057 +0.0026 
ICU.LMO -0.066 -0.12 


0.1, 0.55 


1.57 ±0.14 t° f 9 T 9 


1.5, 4.5 


1 42 ±0 16 + ' 1 +0 ' n 

l.^Z 3ZU.1U -0.096 -0 


4.5, 10.5 


0.92 ±0.08 ±° f 5 t° 078 


10.5, 49 


0.66 ±0.039 + ™ou ^ 


49, 120 


n qc; + n n^7 +0.079 +0.014 

U.OO ZSZV.VOI -0.018 -0.0021 


120, 2000 


0.44 ±0.063 t° l 011 t° T 


85 < E l T < 150 GeV 2 


0, 1 


1.41 ±0.048 1°;? 16 t°oml 


0.1, 0.55 


1.09 ±0.1 ±° ™ T 3 


1.5, 4.5 


n 09 -un 1 +0.095 +0.073 
U.»Z HZU.l -0.032 -0 


4.5, 10.5 


0.68 ±0.057 Toll -o 053 


10.5, 49 


n c;i -un 097 +0.044 +0.033 
u.01 znu.vz,! _o.023 -0.0036 


49, 120 


0.43 ±0.038 j£8» +_ °oll 


120, 2000 


0.41 ±0.048 Toll, Toll, 


150 < E Z T < 700 GeV 2 


0, 1 


0.78 ±0.032 t° T 2 t° T 


0.1, 0.55 


n 79 4-n 1 +0.086 +0.018 
V.IZ 3ZU.1 -0.049 -0 


1.5, 4.5 


0.56 ±0.088 t° °oll ±o° 42 


4.5, 10.5 


n fi7 +Q 07 +0.03 +0.043 
U.UI ZSZV.VI -0.059 -0.028 


10.5, 49 


34 ±0 026 +0 - 015 +0 - 0092 

U.O^ HZU.UZU -0.0057 -0.0098 


49, 120 


n +n n^Q + - 046 + - 016 

U.O^ ICU.UOJ -0.049 -0.0068 


120, 2000 


n 9s -un nAfi +0.019 +0.014 
0.28 ±0.046 _o.oo52 -0.011 



Table 8: Measured ratio R = a(x° hs < 0.75)/V(x° bs > 0.75) as a function of Q 2 

2 

in different regions of E T . 
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Figure 1: (a) Measured dijet cross sections da/dQ 2 /or i" bs > 0.75 (upwards 
triangles) da/dQ 2 /or x° bs < 0.75 (downwards triangles) and da/dQ 2 for the entire 
x obs re gj j0n (black dots). The inner vertical bars represent the statistical uncertain- 
ties of the data, and the outer bars show the statistical and systematic uncertainties 
added in quadrature, except for that associated with the uncertainty in the absolute 
energy scale of the jets (shaded band). The NLO QCD calculations o/DisasterH — h 
(/j 2 = Q 2 + (E^) 2 ) and of Frixione and Ridolfi (fi 2 — (E 3 ^) 2 ) for the photopro- 
duction region are shown for each of the cross sections, (b ) Relative difference of 
the measured dijet cross section da/dQ 2 to the DisasterH — h (fj 2 = Q 2 + (E 3 ^) 2 ) 
and of Frixione and Ridolfi (fi 2 = (E^ 1 ) 2 ) calculations. The hatched band shows 
the theoretical uncertainty of the calculations (see text). 
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Figure 2: (a) Measured dijet cross sections da/dQ 2 for x° bs > 0.75 (upwards 
triangles) da/dQ 2 for x° bs < 0.75 (downwards triangles) and da/dQ 2 for the entire 
x° bs region (black dots). The inner vertical bars represent the statistical uncertain- 
ties of the data, and the outer bars show the statistical and systematic uncertainties 
added in quadrature, except for that associated with the uncertainty in the absolute 
energy scale of the jets (shaded band). The NLO QCD calculations o/DisasterH — h 
with fi 2 = Q 2 are shown for each of the cross sections, (b) Relative difference of the 
measured dijet cross section da/dQ 2 to the DisasterH — (- calculation with fi 2 = Q 2 . 
The hatched band shows the theoretical uncertainty of the calculation. 
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8 9 10 20 30 40 8 9 10 20 30 40 

E J T et1 (GeV) 

Figure 3: Measured dijet cross section d 2 a/dQ 2 dE^ tl (dots). Also shown are the 
NLO QCD calculations of DlSASTERH — h with /j 2 = Q 2 + (E^) 2 and /j 2 = Q 2 , and 
those of Frixione and Ridolfi for the photoproduction region. The scale uncertainty 
for the NLO calculation with n 2 = Q 2 is not shown. Other details are as in Fig. 1. 
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Figure 4: Measured dijet cross section da/dQ 2 dr] F (black dots). The NLO 
QCD calculations of Disaster++ with /i 2 = Q 2 + (E^) 2 and /i 2 = Q 2 as well 
as Frixione and Ridolfi for the photoproduction region are also shown. The scale 
uncertainty for the NLO calculation with fi 2 = Q 2 is not shown. Other details are 
as in Fig. 1. 
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Figure 5: Measured ratio R = a(x° bs < 0.75)/a(a;° bs > 0.75) as a function of 

Q 2 in different regions of E T (black dots). The LO+PS calculations of Herwig 
using the SaS2D photon PDFs are also shown. Other details are as in the caption 
to Fig. 1. 
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Figure 6: Measured ratio R = a(x° bs < 0.75)/a(x° bs > 0.75) as a function of 

Q 2 in different regions of E T (black dots). The NLO QCD calculations of Disas- 
terH — h with fi 2 = Q 2 + (E^ 1 ) 2 and ff 2 = Q 2 as well as the Frixione and Ridolfi 
predictions for the photoproduction region are also shown. The hatched bands rep- 
resent the theoretical uncertainties. Other details are as in the caption to Fig. 1. 
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